Derealization transitions of semi-flexible manifolds 
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Semi-flexible manifolds such as fluid membranes or semi-flexible polymers undergo derealization 
transitions if they are subject to attractive interactions. We study manifolds with short-ranged 
interactions by field-theoretic methods based on the operator product expansion of local interaction 
fields. We apply this approach to manifolds in a random potential. Randomness is always relevant 
for fluid membranes, while for semi-flexible polymers there is a first order transition to the strong 
coupling regime at a finite temperature. 
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Low dimensional manifolds play an important role in 
a variety of different contexts, e.g., as soft matter ob- 
jects or as domain boundaries in condensed matter sys- 
tems. They can perform large shape fluctuations driven 
by entropy Q. According to their fluctuations they can 
be divided into two classes. Flexible manifolds, such as 
interfaces, polymerized membranes, and long polymers, 
fluctuate under a tension controlling their area or length. 
The other class is governed by bending energy; i.e., re- 
gions of high curvature are penalized. Examples are poly- 
mers not much longer than their persistence length, like 
actin or DNA, and fluid membranes. These objects are 
stiffer, and we call them semi-flexible manifolds. 

Whenever a fluctuating manifold is attracted towards 
some other "defect" manifold, there is a competition 
between freely fluctuating configurations favored by en- 
tropy and configurations bound to the defect, which are 
preferred by energy. This competition can lead to a phase 
transition, the so called delocalization or unbinding tran- 
sition. It is often of second order, that is, the amplitude 
of the fluctuations diverges continuously as the transition 
point is approached from within the bound phase. This 
leads to a scaling regime close to the transition whose 
universal characteristics can be described by a contin- 
uum field theory. Well-known examples of delocalization 
are wetting phenomena H. For interfaces and polymers, 
these transitions have been widely studied 0. In the 
case of polymers, even the generalized problem of N mu- 
tually attracting objects can be treated. Analytically 
continued to N = 0, this describes a directed polymer in 
a random medium J4J, which in turn is related to theo- 
ries of stochastic surface growth ||. The delocalization 
transition then corresponds to a roughening transition 
between a smooth and a rough growth mode. 

An important class of low-dimensional manifolds are 
stretched objects with mainly transversal fluctuations. 
These are described by a d-dimensional displacement 
field r(i) which depends on a D-dimensional internal vari- 
able t. The continuum Hamiltonian takes the form 



(1) 



where (V fe r) 2 = £f =1 Y,i =1 {d k n{t)/dt k a ) 2 is the lead- 
ing tension (k = 1) or curvature (k — 2) energy in a 
small-gradient expansion. The potential V describes the 
interaction of the manifold with an external object or 
boundary at r = 0, or the mutual interaction between 
two manifolds with relative displacement r(t). Phys- 
ical realizations of flexible manifolds (k = 1) are di- 
rected polymers and flux lines in a type II supercon- 
ductor (D — l,d — 2) ^,0, steps on a tilted crystal 
surface (D = 1, d = 1), and domain walls in a ferromag- 
net (D = 2, d = 1). Short-ranged interactions of definite 
sign can be represented as V(r(t)) = g§(t) in terms of 
the local contact field 



m ee S(r(t)) . 



(2) 



The scaling dimension of this field, x$ = d\, is given in 
terms of the roughness exponent Xi with \ — (2 — D)/2 
for k = \. There exists a well known perturbative frame- 
work HH to treat such interactions, which has been ap- 
plied extensively to the problem of self-avoiding mani- 
folds Jl0|-[l4[. More complicated short-ranged jl5| and 
long-ranged |16| interactions have been studied as well. 

In this Letter, we develop the field theory of semi- 
flexible manifolds with local interactions, described by 
the Hamiltonian (pi) with k = 2. Physically interesting 
cases are again polymers (D = 1, d = 1,2) fll7|-|l9[ ] and, in 
particular, fluid membranes {D = 2, d = 1) [|20|]. Since a 
semi-flexible manifold has a locally well-defined orienta- 
tion, we have to consider interactions V(r(t), Vr(i)) that 
depend both on the displacement and on the orientation. 
We find there are now two important scaling fields: Local 
contacts at arbitrary orientation are still represented by 
the field <&(£) given by (|j), which has dimension ir$ = d\ 
with x = (4— D)/2. Due to the stiffness, however, closeby 
contacts always have a preferred orientation parallel to 
the defect. Such contacts are described by the field 



n(t) = 6{r(t))S(Vr(t)) 



(3) 




(b) 




FIG. 1. (a) A flexible and (b) a semiflexible polymer r(t) 
bound to an attractive defect indicated by the dashed lines. 
The unbound segments join the defect at arbitrary orientation 
and at fixed orientation dr/dt — 0, respectively. 



ergy density / = F/T D becomes invariant under trans- 
lations of t. This series takes the easiest form if the man- 
ifold is subject to wall constraints forcing the probability 
density p(r') = (S(r(t) — r')} to vanish on the bound- 
ary of the hypercube, in particular along the "edge" 
r = 0. The density then takes the asymptotic scaling 
form p(r) ~ ( ri . . . r d ) e R^ d(1+e ^ for |r| < R, with an ex- 
ponent 9 > expressing long-ranged suppression of the 
configurations close to the boundary. The wall constraint 
is natural in d = 1 for a fluid membrane at a planar sys- 
tem boundary, or for a pair of membranes without mutual 
intersections. The generalization to arbitrary d has been 
chosen such that the Hamiltonian (|lj) remains factoriz- 
able. Short-ranged interactions with the manifold r = 
are now described by the local field 



tt(t) = lim r 

r— ^0 



-dt) 



Ip( r "'(*)- r )> 



(4) 



of dimension xq = d\ + dD{\ ~ !)• The scaling fields 
$ and fi are found to obey an operator product ex- 
pansion. This is a well-known concept in field theory 
(see, e.g., Ref. |yj), which has been applied extensively 
to flexible manifolds p2[. It allows us to write down 
renormalization group equations for generic local inter- 
actions. This leads to results for the derealization of 
semi-flexible polymers and of fluid membranes. Most 
importantly, the bound state of a semi-flexible manifold 
turns out to be maintained by contact interactions (0) 
at fixed orientation, while the bound state of a flexible 
manifold involves interactions of the form (0). Typical 
bound state configurations are compared in Fig. 1 for 
the case of polymers (D — 1). Our one- loop results are 
in agreement with previous results obtained by approxi- 
mate renormalization methods |23| , |24| and by approaches 
specific to polymers ||17Hl9|- They are exact for D = 1 
and can be improved systematically for higher values of 
D. Furthermore, they can be applied to a semi-flexible 
manifold in a quenched random potential (taken to be 
Gauss distributed with mean V(r, t) = and variance 
V(r,t)V(r',t') = a 2 5{t - t')S(r - r')). In the replica for- 
malism, this is equivalent to N interacting semi-flexible 
manifolds in the limit of vanishing N. For fluid mem- 
branes, any amount of disorder is relevant and leads to a 
strong coupling phase. For semi- flexible polymers, how- 
ever, small amounts of disorder are irrelevant (unlike for 
their flexible counterparts). There is now a first-order 
transition to the strong coupling phase at a finite amount 
of disorder. A quantitative description of the disordered 
strong coupling phase is, however, beyond the means of 
perturbation theory. 

We study the manifold displacement field in a hyper- 
cube of longitudinal extension < t a < T (a = 1, . . . , D) 
and of transversal extension < r^ < R (i = 1, . . . , d). 
For T ^S> -R 1 /*, the free energy becomes extensive, F ~ 
T D R~ D / X , and the perturbation series for the free en- 



whose expectation values are finite. Due to the con- 
straint, these interactions are always at fixed orientation 
Vr = 0. Hence, we have used the same symbol Vt as for 
the field (||) of the unconstrained system. The correlation 
functions 

<n(t)) - R- xn/x (5) 

(n{t)n(t'))~\t-t'\- x "(n(t)) + ... (\t-t'\ x ^R) (6) 

define the scaling dimension xq. In the constrained sys- 
tem, the density (0) is linked to the pressure of the system 
by a wall theorem, 

(O) - (df/0R) d - R-dH+D/x) _ (7) 

This determines the exponent value [|25| 

xa = d( X + D), (8) 

which differs from that of the unconstrained system. It 
is in agreement with a conjecture from functional renor- 
malization |g6[ for general D and a direct calculation |27| ] 
for D = 1. 

The interaction part 5f(h, R) = f(h, R) - /(0, R) of 
the free energy of the system with V(r(t), Vr(i)) = hil(t) 
can be expanded as a power series 

h 2 r 

5/ = h(n) + — / d D t(n(o)n(t)) c + o(h 3 ) (9) 

containing the connected correlation functions 

(O(0)fi(f)) c = (fi(O)Q(t)) - (n) 2 (10) 

etc., taken at h = 0. By (0), there is a whole line in the 
{D,d) plane, where the interaction Q is marginal, i.e., 
where xq, — D. This line is given by 



d*{D) = 



2D 
A + D' 



(11) 



The perturbation series (||) has poles in 

e = D - xa = D - d{2 + D/2) , (12) 

which can be regularized around any point on the line 
of marginality ([ll]) ; see the discussion in jl(J for flexible 
manifolds. The singularity of the two-point function (flC 
is given by the operator product expansion 



n(t)fi(t')~ \t-t'\- xn fl(t) 



(13) 



This singularity determines in a standard way |22f| the 
one-loop renormalization group equation of the dimen- 
sionless coupling constant v = hR e ' x . In an appropriate 
scheme, this takes the form 



v = ev-v 2 + 0(v 3 ) 



(14) 



The unstable fixed point v* = e + 0(e 2 ) represents the 
transition. The linearized form v = e*(v — v*) + ... 
with e* = — e + 0(e 2 ) then determines the scaling of the 
transversal localization length £ = (r 2 ) 1 / 2 , 



£~(v*-w)- x / £: 
and the scaling dimension 



(v < v*) , 



D - e* = 2D - x n + 0(e 2 ) , 



(15) 



(16) 



which takes the place of xn in the correlations (H) and 
(H) at the transition point. These relations describe the 
scaling of a bound state maintained by contact forces at 
fixed orientation. Typical configurations look similar to 
those of Fig. 1(b) but are confined to the region r* > 0. 

The most interesting application of (Eq) and (M) is the 
derealization transition of a fluid membrane from a hard 
wall (D = 2, d = 1), where £ ~ (T c - T)~ x (since the ef- 
fective coupling is temperature-dependent) and Xq = 1. 
Not surprisingly, these one-loop results are in agreement 
with those from functional renormalization |23J. They 
also fit very well the numerical values of |24| . which im- 
plies that higher order corrections must be small. The 
system with wall constraint at h = can be regarded 
as an unconstrained system in the limit of large repul- 
sive interaction. Conversely, the scaling at the transition 
point of the constrained system may be related to that 
of the free unconstrained system. Indeed, the one-loop 
value Xq from (um equals the dimension xn — 1 of the 
free field (Q), indicating that the sum of the higher order 
corrections in ( J14] ) and ( fig ) may vanish altogether at the 
specific point (D = 2, d = 1). 

For the case of polymers (D = 1) it is easy to show 
that the multi-point correlations entering (ph factorize 
after "time-ordering" the interaction points, 

<n(ti) . . . n(t n )) = (n(ti))ft(t 2 - *i) • ..n{t n -t n -x) 

(17) 



for ii < . . . < t n . The factors ft(t) = (ft(0)ft(i))/(ft) 
can be interpreted as "return" probabilities. Using ( |17|) , 
one can show that the polymer perturbation series (|9|) 
is one-loop renormalizable; i.e., the connected two-point 
function (Q(0)Q(t)) c generates the only primitive singu- 
larity, and there are no higher-order terms in (|14|). The 
arguments are completely analogous to those for the case 
k = 1 g|||. The implications of @ and @ on general 
semi-flexible polymers are discussed and verified numer- 
ically in p8j . Here we have given a unified derivation of 
these relations, stressing the theoretical analogies with 
their known counterparts for k 



1 |29 2 



We now turn to systems without the wall constraint 
and restrict ourselves to D = 1, namely mutually inter- 
acting semiflexible polymers and, in particular, a single 
such polymer in a random medium. The latter system 
has a possible biostatistical application in the theory of 
sequence alignment ]3(J. In the absence of a wall con- 
straint, we have to study generic contact interactions 
V{t) = g<S>{t) + hn(t) involving the fields (§) and (|). 
The perturbation scries then contains connected corre- 
lations ($(ti) ■ • • $(t„)fi(*i) . . . 0(0) in the free theory 
(g = h = 0), which can be computed exactly. They can 
be shown to obey the operator product expansion 



*(*)*(*') ~ \t-t'\- d n(t) + 

$(t)0(t') 



t-t' 



*0(t) 



\t-t 



r\-2d 



n(t) 



(18) 
(19) 
(20) 



These relations just say that any pair of closeby con- 
tacts amounts to a single contact at fixed orientation, 
multiplied by a singular prefactor. These singularities 
determine the one-loop renormalization group equations 



u = (1 - 3d/2)u 

..2 „,2 



(21) 
v = ev — V — u~ — cuv (22) 

for the dimensionless couplings u = git , ,( 1 ~ 3d / 2 )/x and 
v = hR e / x with 



xn 



1 -2d. 



(23) 



The corresponding flow diagram for d = 1 is shown in 
Fig. pi The unique derealization fixed point (u* = 
0, v* — e) is on the line u = 0. This property ensures 
that the constant c in (E3) drops out of the critical expo- 
nents. It will be preserved for any d > 2/3 also at higher 
orders since any operator product ]1- ^(U)^(t'j) couples 
only to fi(ii). The perturbation scries at u = 0, however, 
is factorizable according to (|l^) and one-loop renormaliz- 
able in exactly the same way as with the wall constraint. 
Hence, the (in D = 1) exact relations (|15| ) and ([16]) still 
hold (with e given by (|23J) and xn = 2d), resulting in 
£ ~ (T c - T) 3 /( 2 - 4 « for 2/3 < d < 1 and x^ =2 -2d. 
This scaling dimension turning negative for d > 1 indi- 
cates that the transition becomes of first order; see the 
discussion and extensive numerics in |p8[ . An analogous 
first-order regime is known for flexible polymers fl29|]. 




FIG. 2. Renormalization group flow of contact interac- 
tions for a semi-flexible polymer in 1 + 1 dimensions. The 
flow equations (Ell), (123) have the single unstable fixed point 
(u, v) — ( — 1,0) marking the derealization transition. 

The above arguments can be generalized to the replica 
theory of N semifiexible polymers coupled by pair con- 
tact forces. If these are of type fl, the time-ordered per- 
turbation series can be mapped term by term onto the 
perturbation series of flexible polymers with $ interac- 
tions. Its leading divergencies are due to "ladder" dia- 
grams with the same pair of polymers interacting at sub- 
sequent points ti |3l]j32| ■ The presence of $ interactions 
for semifiexible polymers does not change these singular- 
ities by the same argument as for N = 2: Due to their 
stiffness, any two semi-flexible polymers interacting twice 
in a short interval have to be parallel to each other, or, in 
other words, the leading divergent diagrams behave like 
diagrams involving only Q operators. For the Q system, 
however, the results of ph,p2| immediately carry over 
and imply that the critical behavior at the derealization 
transition does not depend on N. In particular, the ran- 
dom limit of vanishing N becomes trivial. We conclude 
that a 1 + d dimensional semi-flexible polymer in a ran- 
dom potential has a phase transition between a weak and 
a strong coupling phase at a critical strength of the ran- 
domness for d > 2/3. This phase transition corresponds 
to the roughening transition of the Kardar-Parisi-Zhang 
equation in Ad dimensions. For fluid nrembranes, on the 
other hand, an arbitrarily small amount of disorder is 
relevant and leads to a strong coupling phase. 
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